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ABSTRACT 


The operation of the BBA Real-Time Combiner (RTC) is discussed and its 
performance investigated in detail* It is shown that each channel of the RTC 
can be modelled by Q simple block diagram in the z— transform domain from which 
all pertinent transient and steady state behavioral characteristics can be 
determined. In particular, the characteristic equation of the tracking loop 
and its equiv4^1ent noise bandwidth are found and used to evaluate the closed 
loop transient response and steady-state mean squared timing jitter. The 
impact of the totality of these loop jitter contributions on the combiner 
output SNR is evaluated and illustrated numerically. These results show that 
for parameters of interest to various space missions, the RTC is capable of 
providing significant SNR improvement relative to a single receiving antenna. 
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SECTION I 


INTRODUCTION 


The Deep Space Network (DSN) is responsible for providing interplanetary 
communications for the National Aeronautics and Space Administration (NASA). 

A key feature of the DSN system is the Improvement in effective antenna aper- 
ture gained by combining the outputs of a number of antennae in the network, 
i.e., multiple antenna arraying (see Figure 1). In particular, if system per- 
formance is measured in terms of the ratio of system gain to system noise tem- 
perature (G/T), then the effective G/T after arraying is theoretically equal to 
the sum of the G/T's corresponding to the Individual antennas being combined. 
More realistically, one must have available a detailed analysis of the array 
combining system accounting for all real losses in order to identify how much 
of the above theoretical limit is practically achievable after signal proces- 
sing. The purpose of this report is to provide such an analysis from which 
estimates of performance degradation from the ideal case can be assessed in 
terms of actual system design parameters. The final analytical results will be 
numerically illustrated with several examples corresponding to actual deep 
space mission designs. 

Figure 2 shows the basic structure of the Baseband Assembly (BBA). Its 
purpose is to combine coherently signals from several antennas and perform sub- 
carrier demodulation and symbol synchronization on the combined signal. The 
BBA consists of three main devices! the Real-Time Combiner (RTC), the Sub- 
carrier Demodulation Assembly (SDA) and the Symbol Synchronization Assembly 
(SSA). Ttie last two assemblies are often designated as a single assembly, the 
Demodulation Synchronization Assembly (DSA). The Inputs to the BBA are analog 
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Figure 1. Multiple Antenna Arraying System 
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baseband telemetry signals troni L receivers. These signals are low-pass fil- 
tered, sampled, and digltlKOd. After that, all the processes performed in the 
BBA pertaining to time alignment, subcarricr demodulation and symbol sychron- 
ization may be looked upon as purely digital operations on unitless numbers. 

In this report, the operation ni the RTG will be discussed. The purpose 
of the RTC is to time-align signals from different antennas in order to in- 
crease the strength of the total combined signal. The basic structure of the 
RTC is Illustrated in Figure 3. After low-pass filtering and analog-to-digital 
(A/D) conversion, the signals from L-1 different antennas are cross-correlated 
with the master or reference signal. These cross-correlations produce error 
signals proportional to the relative time delays between these and the refer- 
ence antenna. The error signals are read by the CPU every T, seconds (dur- 
ation of K data symbols), which produces time delays in the First-In-First-Ouc 
(FIFO) devices in order to compensate for the misalignments. Finally, the 
signals are multiplied by weighting factors and summed in the combiner. The 
combined signal goes to the DSA for further processing. 

The analysis of the RTC's performance will be subdivided into the follow- 
ing categories: 

1. Statistics of the Error Signal 

2. RTC Laplace Transform Block Diagram and Transfer Function 

3. Noise Analysis 

4. Delay Rate Compensation 

5. Optimum Coherent Combining 

6. Degradation of the Summer Output Signal-to-Noise Ratio (SNR) Due to 
Delay Jitter 


4 



SIGNAL LPF 11 : J 



5 


Figure 3. Real-Time Combiner 













SECTION r.l 


STATISTICS OF THE ERROR SIGNAL 

Uaseband signals from L antennas are somplect at the Nyqulst rate r^ 
samples/soc. Each sample is converted to an 8-bit binary word by on A/D con- 
verter. The most significant bit of each sample represents the sign of its 
voltage: a "O'* corresponds to a positive voltoge and o "1" corresponds to a 
negative voltage. The A/D output corresponding to the first antenna (herein 
referred to as the master or reference antenna) is cross-correlated with the 

outputs of each of the other L-1 A/D's to form a set of error signals propor- 

th 

tional to the difference in time delay between the reference and the .1 
antenna, j ° 2, 3, ...» L. 

The basic structure of the L-1 identical quadrature correlators is shown 
in Figure Let y^j represent the i*"*^ sample of the antenna A/D 
output and, in particular, y^^ represents the i**'^ sample of the reference 
antenna A/D output. Then, in the upper arm of the correlator, the sign 
bit of y^^ is exclusive OKed (XOR) with a delayed (N^ samples) version of 
the sign bit of y^j^ , while in the lower arm, the sign bit of y^j^ is XORed 
with a delayed version of the sign bit of y^^ . The truth table for the XOR 
operation with arbitrary inputs A and B is given below; 

A B ^ 

0 0 0 

0 1 1 

I 0 1 

1 1 0 


b 










At til output of the XOR gates, the signals and x^^j nssumc values 
•'0" and "I" according to the above table. The sample x,j. Is subtracted from 

L.1,1 

to produce Xg^j, wluch, after proscaling, goes to a counter. The out- 
put oi the counter is coiled the error signal ond Is o random variable whose 
mean and variance are vital to a determination of the RTC's performance. In 

the ensuing analysis, we shall derive these two moments, 
t h t li 

Tne i sample of the j signal, y^j , is modelled as the random 
variable 


” ^ij "ij 


( 1 ) 


where and respectively ^ represent the signal and noise componenta 
of the Input signal# The first and second moments of are given by 
and 

■= ^ 

E(y,/) = E(„y2, + ^ 


( 2 ) 


where represents the amplitude of the signal and o^j is the 
noise amplitude in the channel. Thus, 


rms 


0 . •= N-.B . 

nj Oj nj 


(3) 


where Nqj is the one-sided noise spectral density and is the one-sided 

baseband bandwidth at the BBA input. 

If the signals coming from the reference and the j antenna were noise- 

free and perfectly aligned, then the output of the correlator, x_ . . , would be 

Sij 

zero all of the time. However, noise introduces randomness into the values of 
y^j^ and y^^ , The probability that y^^^ will flip from "0" to "1" or vice 
versa is 
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1 ( r^i \ 

4 crfc [\/ s- 1 

^ \V2a // 

\ ni / 


and the probability that will do the same is 


Pj “ 2 


where erfc x = 1 - erf x is the complementary error function defined by 


r 2 f 

f c X » f 

YV J- 


exp(-y )dy 


Then, in terms of Pj^ and p^ , the probability that or will 


In error is 


PoJ " + '’■-Pj’f'l 


I : - .. 


or, alternately, the probability that Xy^^ or Xy^^^ is correct is 


^■Pej “2 


Given that y^ and y^ are both equally likely to be a '0" or "1", then the 
expected value of xjjij or x^^j is 

= 0 X (1-Pgj) + 1 X Pgj “ Pej’ “ “ y L 


Define s. as follows: 
j 


s. - E « •, <“ 

Then, since random variables a=Uor L, j = 1, 2, ••«, L corre- 

spond to Nyquist samples, they are independent (with respect to the index i) 


and have first and second moments 
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E{8j) « 


( 11 ) 


E{sj^} « NPgj + N(N-l)p^j^ 

and variance 

0 ^ « " (E{sj})^ = Np^j(l-p^j) 

SubsCitutlng Eq. (7) into Eq. (13) gives 



When the above steps are repeated with and y^j having opposite signs, 
the same variance is obtained. Hence, the variance is Independent of the 
transitions of the signal and is the same for deterministic and random data. 

After this preliminary discussion, the means and variances of the error 
signal will be derived for two classes of data: 

1. Random data 

2. Random data plus subcarrier 

A. ERROR SIGNAL STATISTICS FOR RANDOM DATA 

Let Ng be the number of Nyquist samples per symbol and p^ be the prob- 
ability of symbol transition. Because of the delays in the arrival of the 
signal to different antennas and the delays in the correlator, the four streams 
of samples entering the two XOR gates of Figure 4 will have relative displace- 
ments as shown in Figure 5. 

Referring to Figure 5, the expected number of "l"s, in the upper arm aver- 


( 12 ) 


(13) 
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Figure 5. Timing Diagram 
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nged over one symbol time, will be* 


” ^W^e '*■ Pt^U^^-Po^ 


( 15 ) 


where Ny ° N^j “ ° •^he relative delay In the upper arm. Similarly, in 

the lower arm the expected number of "l"s will be 

where ° + 1^ = the relative delay in the lower arm. Subtracting 

Xj^ from Xy , we obtain the error signal, averaged over one symbol, namely. 


^S “ 2Nj^p,j,(2p^-l) 


( 17 ) 


ti h 

Here Nj^ is the relative delay of the two signals during the k symbol. 

Since the number of symbols, K, in an update interval of T seconds is very 

M 

large, via can model the delay error N|^ at the instant k, by a continuous 
variable N(t) (see Figure 6), in which case the average delay error over one 
update interval becomes 


^t+T 
L -^t-T 


R 

N(t)dt 


( 18 ) 


where T. = T + T„. Here T denotes the computation time, i.e., the 

L C K C 

time required by the CPU to process the error signal sample, and Tj^ denotes 
the remainder of the update time interval. 

th 

The expected value of the error signal at the end of one (e.g., the n ) 
update interval will be 


*For simplicity of notation, we shall herein drop the subscripts i and j pre- 
viously used to denote the particular sample arid channel of interest, respec- 
tively. Whenever this is unclear, we shall reintroduce these dependencies. 
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Figure 6. Delay Error Versus Time 
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(19) 


“ I: 


where 1/m^ is the scaling factor at the input to the up/down counter (sec 

Figure 4). Note thot ■=> o when p<p » 0, i.Ci, when no symbol transi- 
tions occur. 

B. ERROR SIGNAL STATISTICS FOR DATA PLUS SUBCARRIER 

Referring to Figure 7, when random date modulates a subcarrier of N 

SC 

Nyquist samples per subcarrier period, then, in the upper arm, the expected 
number of "l"s averaged over one symbol period will be 

+ P.jNuP^ + (l-Pj) for |H„| (20) 


A similar equation is obtained for the output of the lower arm, x^, with Nl 
substituted for N^. 

Subtracting again from xy we get 

where 



( 22 ) 


Note, that for large values of N /N (the data period is much larger than 

s s c 

the subcarrier period), the error signal is Insensitive to the probability of 
symbol transition. This is true in the region jNj^j £ ^sc'^^ only. The 
error signal at the end of one update Interval will be 


; = ICN (2p ~l)f 7ra 
n e ^e p o 


( 23 ) 
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Comparing Eq. (19) with Eq. (23), we can express the expected volue of the 
error signal ot the end of one updote interval for both the random data nnd 
rondom data v;lth subcarricr cases os 


e ^ G ’ N 
n Q e 


(2A) 


where is the gain which relates the overage error signol at the end of an 
update interval to the overage delay error over thot Interval ond is defined 

by 


« -J^erf 
"’o 


erf 




2a 


la 


2 \ ^p 


(25) 


'nl / \ ■ -nj 

with fp “ 2p^ for random data without subcarrier and jN^j N^, and 


f = f’ = 2 
P P 


( 2N \ 

N~^ “ p,j,l for random data plus subcarrier and |N^ | ^ Ng^/2, 


Eq. (2A) can be rewritten (using Eq. (18)) as 


where 


a 


n 


G y N(t)dt 


(26) 



(27) 


Herein, 


will be Created as the 


"gain" of the cross-correlator-counter. 


The variance of the error signal does not depend on the symbol transition 
rate and equals the sum of the variances in the upper and lower arms. Using 
Eq. (13), the variance at the output of the non-correlator counter will be 
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The above equation applies to random data, random data plus subcarrier, or a 
squarewave* 
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SECTION III 


RTC LAPLACE TRANSFORM BLOCK DIAGRAM AND TRANSFER FUNCTION 


The throe basic components of the Real-Time Combiner are*, cross-correla- 
tor-counter, loop filter and FIFO shifter. The transfer functions of these 
components are modelled as follows:*^ 

CrosB-Correlator-Counter 


Loop Filter 


FIFO 


e(z) A 
Ng(s) “ 


C(z.s) = Gq i 


z 


AN (2) A 
e(z) " 


L(z) 




zA 

2-B 


N(s) A 
AN(2) “ 


F(z,s) 



S 2-1 


<29a) 


(29b) 


(29c) 


The time delay between the time instants that the error signal is 
read by the CPU and the delay update AN^^ produced is modelled as 


-gT^s 

D(s) » e 


(30) 


where g = Ic/T^ is the ratio of computation time to loop update time. Al- 
though D(s) is between the cross-correlator-counter and the loop filter, in 
order to obtain the equivalent block diagram in the z-domaln we place D(s) 
after the FIFO. 


*The transfer functions are expressed as either z-transforms or hybrid s/z- 
transforms as appropriate, depending on whether the operation being performed 
is purely digital or a combination of analog and digital. 
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Combining ISqs. (27) through (30) • the mixed mode Laplace block diagram oi: 
Figure 8 is obtained. The parameter C in this figure is the gain associated 
with delay rate compensotion. Referring to Figure 8 


X(s) 


N(s) 

8 


Y(z) 





(31) 


T,^R 

But since z corresponds to an update Interval of T)^ seconds, i.c,, z « e 
then 


1 




(32) 


Thus 




X(s) 


8 


-ST. 
e ^ AN(z) 


(33) 


v/here AN(z) is the z transform of the FIFO input. Taking the z~transform of 


Eq. (33), we obtain 



whore the asterisk denotes the z-transform operation. 


(3A) 


In general 

r >gT 8 y 

e R(s) A R(z,m), m “ 1-g (35) 

where R(z,m) is the modified z-transform of R(s). Using this property, 

Eq. (3A) becomes 

X(z) * - ^N(z) (36) 

With this transformation, Figure 8 has the equivalent block diagram of Fig- 
ure 9, which is entirely in the z-plane. 
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Figure 9. RTC Equivalent Block Diagram in z Domain 
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TliG open-loop transfer function for tlie RTC la 


Hl<» ) 

X(z 


L A G(z^ o si 

8) o (z-1) (z-B) 


( 37 ) 


where G la the closed loop gain which for this configuration becomes Identical 
with Gp as defined in Eq. (27). The closed-loop transfer function for the 


RTC is 


or 




P 




H(z) 


G(z) 

1+G(z) 


(38) 


H(z) ■= AG -y 


z(l-s) + R 


z^ -!• z(AG(l-g) - B - + B + gAG 

The denominator of Eq. (39) the loop characteristic equation and has the 
root locus diagram as illustrated in Figure 10. Expressing it as 


(39) 


C(z) => z + a^z + 8 q 


(40) 


and solving for the gains A and B, we obtain 


ao + + 1 


(41) 


B = aQ(l-g) - g(aj^+l) (42) 

For given values of damping factor ? and settling time t^, the natural fre- 
quency of the loop will be 



(43) 


The values for the constants a^ and can be found from the following 
relation 
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- NOMINAL GAIN, AG = 0.2374 
#3 - AG = 1.8 


Figure 10. Root Locus Diagram 


23 



o - 2 e 


( 44 ) 


a 


1 


n L 


cos 



-25W T. 
^ n Tj 


(45) 


where ag&ln Tj^ is the loop update time in seconds. This selection of Bq, 
Sj, and, consequently, the gains A and B will Insure that the poles of the 
transfer function H(z) will lie inside the unit circle and the loop will be 
stable. For example, by making ? = 0.707, t = 9.5 sec., and T, = 1.0 
sec., we get Sq ° 0.4308, “ -1.1934, and the natural frequency * 

0.6 rad/sec. Also, the poles of H(z) will be located at 0.5967 + j 0.2734. 
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SECTION IV 


NOISE ANALYSIS 


The one-sided noise-equivalent bandwidth is defined as 




L H (1) 

Letting z = 1 in Eq. (39) gives H^Cl) = 1. Defining the integral I„ as 


( 46 ) 


f 

2iTj / 


•^H(z)H(z'^) , 


(47) 


we obtain 

((1-b)^ H- R^)(l + AGr •{• B) - 2(l-g)g(AGq-g) - B - 1) 

^ (1 - (AGg + B)^) (1 + AGg + B) - (1 - (AGg + B)) (AG(l-g) -B-1)^ 


The one-sided noise-equivalent bandwidth now becomes 

1 

“l 21V 2 

U 


(49) 


Table 1 gives for selected values of AG and g. One can observe from these 
results that B is almost insensitive to the parameter g. 

i./ 

Referring to Figure 9> the open loop noise process n at the output of 

the cross-correlator-counter has zero mean and variance given by Eq. (28). We 

want to find the variance, o^, of the closed loop delay error as a func- 
2 

tion of 

e 

Assuming that the spectrum of the noise process, n, is wide in relation 
to the loop bandwidth, B, , the closed-loop variance of the steady-state error 
signal will be 




H(z)H(z 



(50) 
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Table 1. Noise Equivalent Bandwidth in Hz vs. Loop Gain AG 
and Computation Time Factor g 


Loop Gain 


Computation Time Factor 


AG 

g = 0.0 

g = 0.3 

g « 1.0 

0.0001 

4.4 X 10“5 

3.9 X 10"5 

3.1 X 10-5 

0.0010 

4.4 X 10“^ 

4.0 X lO"'^ 

3.1 X lO"*^ 

0.0100 

4.4 X 10“^ 

3.9 X 10“3 

3.2 X 10“3 

0.1000 

4.6 X 10"2 

4.12 X 10"2 

3.7 X 10"2 

0.2374 

0.1137 

0.1058 

0.1133 

0.5000 

0.2661 

0.2630 

0.4716 

1.0000 

0.6751 

0.7851 

— 
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v^here U(z) is given by Eq. (39). Comparing Eq. (50) with Eq. (46), we see that 
the closed Loop variance, can be expressed as a function of either 

Bl or l 2 » Since H(l) ° l.O, then 


“c ^ , \ (51, 

8S 

2 

The closed loop variance, Oj^, of the delay error is readily obtained from 
the steady-state closed-loop variance of the error signal, namely, 


0 


2 

N 



(52) 


Combining Eqs. (25), (^7), (28) and (51), and substituting the result in Eq. 

(52), we obtain, after much simplification, the desired result, namely the 

t ll 

variance of the delay error in the j RTC, i.e. , 



B-T^N 
L L s 


Kf, 


erf 



(53) 
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SECTION V 


DELAY RATE COMPENSATION 


Referring do Figure 8, we want to find the value for the gain C such that 
when an estimated delay rate compensation Is applied at the input of C, the 
expected value of the error signal Is zero. The estimated delay rate can be 
expressed as 


n^(t) = (n + An) u(t) ( 54 ) 

where u(t) is the unit step function, n Is the actual delay rate of the Input 
signal and An is the error made in estimating n. Neglecting for the moment 
the noise process n, the state equation of the RTC loop is 


( 


«i+l 

AN 


1+1 


) 



-AGg - B 
AGg + B , 



/(Tj^-Ac(l - 2g) - C) 
\(Ac(l - 2g) + C) 


- C 


C 



(55) 


where N^ is the delay error at time 1, AN^ Is the delay update at time 1 
and c = GTl/2. 

In the steady state (1+®), we want 


" “ 2 


(56) 


Inserting Eq. (56) into Eq. (55) and solving for N, we get 


N e - 




(57) 


Inserting now Eq. (57) in the upper or lower row of Eq. (55) and solving for C, 
we get the expression for the doppler compensation gain, namely, 


C = T^ (1 - B) 

Xj 


( 58 ) 





Note that the gain C is Independent of the computation tline> 

Taking the z~transform of N(b)/ 8, the steady state error due to an input 
delay ramp when n(.(t) = 0 is 

/!a^h(.z-n)\ / A o^(z-i) 

I 2 (.- 1)3 ) W 


nTj^(l - B) 


(59) 


Using Eq. (38) when N(s)/s = 0 and taking the z-transform of Eq. (5^), the 
steady state error due to delay rate compensation n^(t) is 

- Hn, ^^A4.AA^ -^1 (-C) G(z(l-g) + fl) 

" U ^ z-lj ^ z ^ (z-l)(z-B)[l + G(z)] 


■bs2 


“ -(n+An) ” (60) 

Substituting Eq. (58) into Eq. (60) and using Eq. (59), we get the total steady 
state error 

%s ~ ^ssl *^ss2 

nTj^d - B) hTj^d - B) AnTj^(l - B) 


AnTj^(l - B) 
A 


(61) 


Since E{An} = 0, the steady state error has zero expected value, as desired. 
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SECTION VI 


OPTIMUM COHERENT COMBINING 


In this section, we determine the optimum weighting factors, a^; 1 ° 1, 
2, ..., L, such that the SNR at the combiner's output is maximized (see Fig- 
ure 3). Referring to Figure 11, coherent signal combining implies that the 
amplitudes of the signals add coherently. This means that at the combiner 
output 



( \ 

( L 1 

L 


A 

{^}= IE. 


= E 

j=l ^ ^ 

(62) 


where s^j = the signal amplitude of the i*^^ sample of the signal ® \/^ 

th 

on the average, = the weighting factor for the j signal, and L = the 
number of signals combined. 

The noise processes in each of the channels are statistically independent. 

th 

As a result, their powers add, i.e. , the noise power of the i sample at the 
combiner's output is equal to the weighted sum of the noise powers at the 
input. Hence, 


N 


A 






(63) 


^ f-h 

where } Is the rms noise amplitude squared in the channel* 

The SNR at the output of the signal combiner for an arbitrary set of 
weighting factors is 


R 




2 : 
a, a . 
j nj 


(6A) 
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Given a BOt of signal amplltudoB, noiso atnpllcudosi 

at the comblnor's input, then (a^} is chosen so as to maximize the parametct 
R of Eq. (64). Differentiating R with respect to each of the weighting factors 
and equating the results to zero, i.c., 


3R 

3a, 


» 0 for j “1, . . . , L 


(65) 


provides the solution 


j 




( 66 ) 


Since the common factor does not affect the combiner output 

SNR, it can be set to any arbitrary value, c.g., unity. This, the equivalent 
choice of weighting coefficients, is 


j 


2 

nj 


(67) 


This optimization assumes that there are no delay errors at the input to 
the combiner, i.e., the BBA's tracking loops are functioning ideally. 
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SECTION VII 


DEGRADATION OF THE SUMMER OUTPUT SIGNAU-TO-NOISE RATIO 
(SNR) DUE TO DELAY JITTER 


For on orblcrary sot of weighting coef detents , {Oj}> J “ 1» 2 , 
Lt the output of the combiner, ns shown in Figure 11, Is given by 


>‘ai “ “l “jl + J?1 “j "jt' ^ 


2 , • • • » Ng ( 68) 


where j » 1 refers to the reference channel whose signal is assumed to have no 
timing error. Adding oil of the samples over the symbol Interval, 
the output of the combiner In 

^/ \ 


“j v/^'*'s ■ ^2 “j'*3 


- - z ('in - + 1) + ces) 

■%> ^ 

where, as before, n = N /N _ = the ratio of the data period to the subcorrler 
period, K = the number of symbols between deloy updates, and » the delay 
error (number of samples) in the channel during the k*"*’ symbol with 
respect to the reference (1®*") channel. By definition, f°>^ oil k “ 1, 2, 

• t » I K • 

The quantity \/i^^ is the noise at the combiner output over the interval 

tl ll 

of one (the k ) symbol. It represents the noise effects from all channels. 

It is modelled as a Gaussian random variable with zero mean and variance 


2 ^2 2., ,, „ 
0. = X «.o , N = H, N 
1< J nj s A £ 
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Furthermore, bj rcprcBontB random data from adjacent symbols, which 
takes on values +1 with equal probability. Thus, 

E |(4n - bj + ijj » 4n + 1 S u 

< 

E |(4n + 1 - bj)^| » (4n)^ 8n + 2 4 w 
It is assumed that each delay error can be modcllea as a Gaussian 


random variable with mean 


«jk “ "j + 


and variance 


^jk 


E I [nj H- - (Sj + Ajk)] 


Here o„ ^ is the variance of the delay error at delay update instants 

and is given by Eq. (53). Also, nj denotes the delay rate of the j ‘ signal 
in units of samples per symbol. 

Since we have assumed that the delay error is a Gaussian random variable, 
the probability density function of the absolute value of the delay error will 


P(l»jkl) 




with mean 


^jk “ “ \[^ 


2 2 
/ 2 -h^ ^ -c \ 

f - - - 2 H" Jjk 


+ (erf (b) - erf (c)) ~ 


34 



where 


" ° ' 


"s - ‘jk 


N +6 

C o -f- y cl,. » ” (erf (b) + erf (c)) 


A 

Here, “S,. and ^ are defined by Eqs. (72) ond (73), respectively. 


2 2 

The expected value of 6.. ^ 


^be***^ + ce"*^ 

Jk « jk 6 1 


,yi„2 

‘jk V /J“jk 


Assuming zero differential dcipler, l.e. , nj ° 0 for all j, 


a B 0 , b = c 


where 




A —J. 


i “ 


and Eqs. (75) and (77) reduce to 


-l/2p. 


Pj^ 
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(81) 


.2 A a 2 ^ „2 


1 - 


j 


Note that as l/pj appcoachos zero, 


-i/2pj 


" v-j ' 


and for 1/pj that is very large, 


2 




.2 2 

j 


(82) 


(83) 


The SNR at the output of the combiner averaged over the interval of one 
delay update period is given by 



To compute the SNR of the combiner's output, only the case of the zero 
differential doppler will be considered. For this case, both the numerator 
and denominator of Eq. (84) are not functions of k. Hence, Eq, (84) sim- 
plifies to 


R, 



N 

J J s 


Z a (4n + 1 

J J J 


N.N 
A s 




(85) 
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Expanding the numerator of the above equation, we obtain 


Evaluating the expectations os required In Eq. (86) gives 

L _ L 

a . A . + w 

j=2 ^ j«2 






( 86 ) 


(87) 


whore A^ and A^ ore defined by Eqs. (80) and (81), respectively, and u and 
w are defined by Eq. (71). Using these results, the SNR at the output of the 
combiner will bo 


S.N 

R = 

A N. 


L A, 2 

Sa m2 
A N 

s 


( L. vty, a ^ o 

2u X) a - w £ «A 

1=9 J v/r "c j=2 


j=2 J /S^ s 


+ „2 V A V „ -5 a 

h i h N? V h ^ ^ "s 


V n2 


s 


S. N 
A s 


(88) 


where 

D A 1 - A (89) 

is a degradation factor representing the reduction in output SNR from its value 
SaNs/Na when there is no delay error. Hence, at zero delay errors, D = 1. The 
variation about D = 1 is given by 
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2 


( 90 ) 


S. A 


, A. 1* n Uj 

A“2u £ a -^^-w o. 

j»2 ^ ”s j“2 ^ “a N^ 


2 V 2 !i / 4; 

u H s “ (u I: a 

j«2 ^ A N_ \ j«2 


/r A 



The above expression is valid Oor an arbitrary set of (“j}* the optimum 
values are substituted from Eq, (66), then A^^^ is 



then the optimum weighting factors can be specified as 


= 2 , 


1 for j = 2, . . . , L 


Under the above assumptions 


(93) 



(9^ia) 


= L + 3 
nj 


(94b) 
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1 

L + 3 


for j « 2, . . . , L 


(9Ac) 


Finally, it can be assumed that the rms delay errors in all tracking loops 


are Identical, l.e., 


and furthermore, 


“ Oj^ for j = 2, L 

J 


Pj = p for j = 2, 3, . . . , L 


Making all of these substitutions In Eq. (91) gives the simplified result 


« /2 L - 1 

AU J— z — ; — r p r 
V r L + 3 * c 


L - 1 / fT 


22 L-1 22 22/L-l 


(L + 3) 


7T P r - u — - — r-T P t 
2* (j 7r(L + 3 o 


where 




-l/2p^ 


erf(^) 

^p/2^ 


The SNR degradation from ideal performance is given by 

D . = 1 - A ^ (99! 

opt opt 

Using Eqs. (97) and (99), and the above assumptions, the SNR degradation 
at the combiner output is plotted in Figures 12 through 17 versus the input 
symbol SNR of the reference channel for parameter values typical of the Voy- 


ager, Pioneer, and Venus Radar Mapper (VRM) missions. 
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DEGRADATION FACTOR 



REFERENCE ANTENNA SYMBOL SNR, dB 


Figure 12. Degradation Factor vs. Symbol SNR, Voyager Mission, L = 2. 

Other conditions include: noise bandwidth (B„) = 3.75 x 

10^ Hz, symbol rate (r) = 40 x 10^ symbols per second (sps), 
subcarrier frequency = 360 x 10^ Hz, integration time = 
i second, g = 0.3, and = Sj/4; i = 2 L. 





DEGRADATION FACTOR 



Figure 13. Degradation Factor vs# Symbol SNR, Voyager Mission, L = 4. 

Other conditions include: noise bandwidth (B^) » 3*75 x 

10^ Hz, symbol rate (r) = 40 x 10^ sps, subcarrier 
frequency =» 360 x 10^ Hz, integration time « 1 second, g 
0.3, and Sj[^ ^ S|/4^ i ~ 2, L# 
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DEGRADATION FACTOR 



REFERENCE ANTENNA SYMBOL SNR, dB 


Figure 14. Degradation Factor vs. Symbol SNR, Pioneer Mission, L = 2. 

Other conditions Include: noise bandwidth (Bj^) = 135 x 

10^ Hz. symbol rate (r) = 8 sps, subcarrier frequency = 

2 X 10^ Hz, integration time = 1 second, g = 0.3, and S^ = 
Sj ^/45 i “ 2, ..., L. 
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DEGRADATION FACTOR 
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REFERENCE ANTENNA SYMBOL SNR, dB 


Figure 15. 


Degradation Factor vs. Symbol SNR, Pioneer Mission, L = 4. 
Other conditions include: noise bandwidth (B^) = 135 x 

10^ Hz, symbol rate (r) = 8 sps, subcarrier frequency ® 

2 X 10^ Hz, Integration time = 1 second, g “ 0.3, and S^ = 
S^/4^ i ^ 2, ..., L. 
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DEGRADATION FACTOR 



REFERENCE ANTENNA SYMBOL SNR, dB 


Figure 16. Degradation Factor vs. Symbol SNR, VRM Mission, L = 2. 

Other conditions include: noise bandwidth (B^) = 3.75 x 

10^ Hz, symbol rate (r) = 550 x 10^ sps, subcarrier frequency 
= 900 X 10^ Hz, Integration time = 1 second, g = 0.3, and Si = 
S|^/4j i “ 2, ..., 1.. 
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degradation factor 



REFERENCE ANTENNA SYMBOL SNR, dB 


Degradation Factor vs. Symbol SNR, Pioneer Mission, L = 4. 
Other conditions include: noise bandwidth (B^j) = 3.75 x 

10^ Hz, symbol rate (r) = 550 x 10^ sps, subcarrier frequency 
= 900 X 10^ Hz, Integration time = 1 second, g = 0.3, and S^ 
S^/4j i “ 2, ..., li. 
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